Interplay between quantum criticality and geometrical frustration 
in Fe3Mo3N with stella quadrangula lattice 
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> ' Abstract 

ov 

In the ?7-carbide-type correlated-electron metal FesMosN, ferromagnetism is abruptly induced 
to it or when it is doped with a slight amount of impurity (~5% Co). We observed a peak in 



from a nonmagnetic non- Fermi- liquid ground state either when a magnetic field (~14 T) applied 



the paramagnetic neutron scattering intensity at finite wave vectors, revealing the presence of the 
antiferromagnetic (AF) correlation hidden in the magnetic measurements. It causes a new type 
of geometrical frustration in the stellla quadrangula lattice of the Fe sublattice. We propose that 
the frustrated AF correlation suppresses the F correlation to its marginal point and is therfore 
responsible for the origin of the ferromagnetic (F) quantum critical behavior in pure FesMosN. 

PACS numbers: 75.30.Kz, 74.40.Kb, 75.10.Jm 



The geometric frustration in electronic degrees of freedom is a key factor responsible for 
the suppression of long-range magnetic order that may lead to the realization of new and 
exotic quantum phenomena l| . Although the role of frustration was established in localized 
electron systems, no consensus has been achieved as far as their role in itinerant electron 
metals is concerned |2J. One of the reasons for this is the fact that a limited number of 
materials of this type are available. Y(Sc)Mn2, /3-Mn, LiV204, etc. are some of the known 
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examples [3|, whose electronic quasiparticle mass is commonly enhanced but the manner 
in which it is enhanced is still under debate. To resolve this issue, it is essential that 
we find other materials of a similar nature. In this letter, we propose the //-carbide-type 
transition metal compound as a good candidate for testing geometric frustration in metallic 
magnets; the 7/-carbide-type structure (prototype: Fe 3 W 3 C, space group: Fd3m) includes 
the stella quadrangula (SQ) lattice (Wyckoff positions 16c? and 32e) [4], which is a corner- 
shared network of stellate (tetra-capped) tetrahedrons (see Fig. [p. In an SQ as a unit, a 
small regular tetrahedron at 32e is nested in each I6d tetrahedron having the same center 
of gravity. Furthermore, direct coupling is expected at the relatively short Fe-Fe bonds. 
Even though 77-carbide-type transition-metal compounds have been known for a long time 
as impurities in steel and refractory and hard materials, as well as promising candidates 
for catalysts, etc. fl, their quantum physica. properties have been .ess studied thus far. 
(as an exception, see [6].) We have recently reported that FesMosN with the ^-carbide- 
type structure [2j exhibits a non-Fermi liquid (NFL) behavior near the ferromagnetic (F) 
quantum critical point (QCP) under ambient pressure [8|, which is different from the findings 
of literatures published previously [9|, Il0j]; C/T (C: specific heat, T: temperature) shows a 
-logT divergence and reaches 128 mJ/(f.u.mol K 2 ) at 0.5 K, and the resistivity shows a T 5 ' 3 



dependence. FesMosN is a one of the ideal Fe-based materials exhibiting NFL behavior 11 ] . 

The susceptibility x °f FesMosN obeys the Curie- Weiss (CW) law at high T (effective 

moment p c ^ = 2.14 /xg/Fe, Weiss temperature 6 = 2 K) and exhibits a broad maximum at 



^75 K 



Because the x{T) maximum is commonly observed among exchange-enhanced 



Pauli paramagnetic metals undergoing itinerant-electron metamagnetism (IEM) [12J], the 
application of a high magnetic field to FesMosN may yield interesting results. IEM of 
FesMosN is also fascinating related with the recent interest in the phase diagram in the 



visinity of F-QCP and tuning of quantum criticality 13]. While we have detected the 



IEM at ~14 T, it is somewhat different from a typical IEM 12] . In order to determine 



the reason for this difference, impurity doping (Co for Fe) and inelastic neutron scattering 
measurements of FesMoaN were performed. In this letter, we discuss the characteristics of 
the SQ lattice as a new geometrically frustrated system, and we use this model to explain 
the absence of long-range order in pure FeaMosN in spite of the strong electron correlation. 
Polycrystalline samples of (Fe 1 _ x Co x )3Mo3N, where < x < 1, were synthesized by a 



solid-state reaction of a mixture of transition metal oxides in a H2-N2 mixed gas stream [14]. 
Fe203, C03O4, and M0O3 were mixed in a molar ratio of (1 — x)/2 : x/3 : 1 and placed in 
a silica tube; and the mixture was then fired in a gas stream of N2 containing 10% H2 at 
700°C for 48 h followed by heat treatment at 1000°C for 48 h. In order to homogenize the 
samples, the heat treatment, interspersed with intermediate grinding, was repeated at least 
four times, resulting in a systematic variation in the magnetism against x, which is not in 



accordance with the literature data lOj. The magnetization M at low fields was measured 
using a SQUID magnetometer, MPMS (Quantum Design) equipped in the LTM center, 
Kyoto University. The high-field magnetization for field strengths up to 54 T was measured 
for pure FesMosN using a pulse magnet equipped in ISSP at 4.2-100 K. Polycrystalline 
powder was filled into a cylindrical polyethylene tube that measured 6 mm in length and 2.5 
mm in diameter. For the neutron scattering experiments, approximately 20 g of the powder 
was packed into a vanadium tube with a diameter of 20 mm. A triple-axis spectrometer 
ISSP-HER installed at the Cl-1 cold neutron guide of the research reactor JRR-3M at Japan 
Atomic Energy Agency (JAEA), Tokai, was employed. All measurements were performed at 
a fixed final wave vector kf = 1.45 A -1 , with a collimation sequence of guide-40'-open-open, 
and with a horizontal focusing analyzer. The energy resolution was estimated to be 0.22 
meV at zero energy transfer. 

Figure [2] shows the results of the high-field magnetization measurements. The signal 
from the pulse magnet corresponding to dM/dH (at 4.2 K) is plotted against the exter- 
nal magnetic field H in the inset of Fig. [2](a). dM/dH shows diverging behavior at ~14 
T, corresponding to a sharp jump in the M-H curve, which was obtained by integrating 
dM/dH. The metamagnetic fields Hq are different for the field- increasing and decreasing 
processes because of the first-order transition, although the difference (hysteresis width) AH 
is much smaller (0.56 T at 4.2 K) than that for a typical IEM [l2J. AH decreases with T 
and vanishes at Tqm — 40 K (Fig. [2(c)), which is the temperature at which the first-order 
IEM disappears, and corresponds to the substantial magnetic interaction in this compound. 



The increase in Hq is roughly proportional to T 2 up to Tcm (Fig. Mjo) closed circles), as 
commonly observed in IEM 12J. In the field-increasing process, a small step was observed 
at a field (defined as He) slightly grater than H c . H E decreases gradually with T and dis- 
appears at ~20 K (Fig. J2Tb) open circles). The magnetization jump is sharp unlike that 
observed in other IEMs 12[ and rather similar to that of a spin flip in a localized spin sys- 
tem. Successive magnetic transitions in the magnetic field are frequently observed in local 
spin systems; however, they are unlikely to be a part of a conventional IEM with uniform 
spin polarization, thus suggesting the presence of competing interactions or multiple order 
parameters. 

The effect of impurity doping was investigated by replacing Co with Fe. We have success- 
fully synthesized a solid solution between FesMoaN and the Pauli paramagnetic C03M03N; 
the lattice parameter was found to vary linearly with the Co fraction x (Vegard's law) 
in accordance with the literature data [lOj|. As a typical example, the magnetic data of 
(Feo.8Coo.2)sMo3N are presented in Figs. [31(a) and (b). At high T, x follows the CW law 
well (peff = 2.55 yU B /Fe, 9 = 17.1 K). The Arrott plot (M 2 against H/M plot) shows good 
linearity and intersects the vertical axis below 22 K, indicating the occurence of spontaneous 
magnetization p s (=0.28 /xb at 5 K) below the Curie temperature T c = 22 K, despite the 
fact that both end compounds have nonmagnetic ground states. The value of p e g/p s is much 
larger than unity, as is characteristic of a typical weak itinerant electron ferromagnet. Simi- 
lar analyses were carried out for other specimens to enable the construction of the magnetic 
phase diagram shown in Fig. 131(c) (detailed magnetic data will be published elsewhere). The 
characteristics of this phase diagram are as follows. A ferromagnetic phase appears by slight 
doping (x < 0.05). With an increase in x, Tq increases rapidly up to a maximum at x ~ 0.2 
and decreases almost linearly after the maximum up to x ~ 0.7. It should be noted that 
the Tcm of pure FeaMoaN is much higher than the maximum value of Tc and is smoothly 
connected to the variation in Tc for x > 0.2. In other words, long-range ferromagnet ism 
appears to be suppressed particularly at x ~ 0, in spite of strong magnetic correlation. 

To obtain information on the magnetic interactions in FesMosN, inelastic neutron scat- 
tering experiments were performed with the powder; unfortunately, single crystals have not 
been obtained. Figure H] shows the wave number (K) dependences of quasi-elastic neutron 
scattering measured at 5.5 K. Relatively strong scattering was observed, centered at finite 
but relatively small K, which roughly corresponds to the mean distances between Fe atoms, 



suggesting the presence of AF spin correlation together with the F correlation that was de- 
tected from macroscopic measurements. The .E-scan spectra, shown in the inset, were fitted 
with a Lorentzian quasi-elastic scattering function. The spectral widths, i.e., characteristic 
energies of spin fluctuations, are obtained as 0.4, 0.3 and 0.9 meV for K = 1.22, 1.78, and 
2.83 in the reciprocal lattice unit, respectively. 

Let us understand the geometrical characteristics of the SQ lattice on the basis of the 
Heisenberg model. Although FesMoaN appears to behave like an itinerant electron magnet, 
we believe that its spin density is relatively localized at the atomic site, as deduced from 
the strong electron correlation. The SQ lattice consists of lQd and 32e sites arranged in 
the Fd3m space group. In actual 77-carbide-type compounds, the nearest-neighbor (nn) and 
next-nearest-neighbor (nnn) distances are close and become exactly the same when z = 0.3, 

n 

where z is the coordinate of the 32e site (z, z, z). Experimentally, z = 0.2937 [7J and similar 



values [15(| have been reported for FesMoaN, where the nn 16<i-32e distance is 2.387 A and 
nnn 32e-32e distance is 2.549 A at room temperature. Here, we consider only these nn and 
nnn interactions, namely, J\ and J 2 interactions (see Fig. [TJ), respectively, and neglect the 
interactions of the neighbours located further along. The signs of J\ and J 2 are likely to vary 
because, empirically, both the F and the AF states appear in doped Fe-based 7/-carbide-type 
compounds [6j. First, let us understand the nature of an isolated SQ. In an extreme case, 
with J\ < and J 2 = 0, the SQ is not frustrated because the AF spins can be alternately 
assigned to the eight atoms. Needless to say, the same is the case for J\ > and J 2 = 0. On 
the other hand, when Ji = and J 2 < 0, the isolated 32e tetrahedron is a typical frustrated 
unit. These facts suggest that the SQ is frustrated only when J 2 is negative and dominant. 
It should be noted that this is true not only when J x < but also when J x > 0. To verify 
whether this is true for an infinite SQ lattice, we calculated the dispersions of the Fourier 
transform of the exchange integral matrix J(q) among 12 Fe atoms in a unit cell by assuming 
various J\j J 2 ratios. Two typical results are shown in Fig. El A flat dispersion is found along 
the highest branch when J\ is not dominant (Fig. ETa)), suggesting the degeneracies of the 



ground state, i.e., the presence of geometric frustration [16]. On the other hand, J(q) takes 
a maximum at the T point, where the magnitude of J\ is sufficiently larger than that of J 2 
(Fig. EJ^b)), indicating that J\ tends to release the frustration. We confirmed that similar 
dispersions were obtained independent of the sign of J\. Thus, the SQ lattice is a unique 
frustrated system in which frustration is controlled by the | Ji|/J 2 ratio but not affected by 



the sign of J\. 

The anomalous magnetism of FesMosN can be easily understood on the basis of these 
characteristics. The experimental results indicate the coexistence of the F and the AF 
interactions, which are formally denoted as J\ and J 2 , respectively, as one of the simplest 
possibilities. The suppression of long-range order in pure FesMosN is due to the frustration in 
the SQ lattice (J2 dominant case). In other words, the frustrated AF interaction suppresses 
the F correlation to its marginal point. Metamagnetism can be interpreted as the transition 
to the J\ dominant state, assisted by the external field. The onset of static magnetism by 
slight doping is one of common features of a frustrated system. It is not surprising that 
the F order is stabilized when the F interaction coexists with the AF interaction. We have 
tentatively calculated the spin correlation function S(K, E) for the SQ lattice, assuming 
(1) J\ > and dominant J2 < 0, (2) 2 up and 2 down spins in the 32e tetrahedron, (3) 
an up spin at a 16d site when the neighboring 32e spins are 2 up and 1 down, and vice 
versa, (4) no correlation between different SQs, and (5) the magnetic form factor of Fe 3+ . 
Some of these assumptions may not be valid for the actual case, but have been made for 
the sake of simplicity. S(K, E) thus calculated and, furthermore, averaged over equal K is 
illustrated by the solid curve in Fig. HI Although S(K, E) should ideally be compared with 
the E averaged values, there is reasonable agreement between the experimental results and 
the calculations carried out using our model, indicating that our model is consistent with 
the experimental results. Although we have applied the localized moment scheme here, we 
have to reconcile the Heisenberg frustration and the electron itinerancy in further studies. 

In summary, we found a steep metamagnetic transition in the 77- carbide-type Fe-based 
compound FesMoaN from the nonmagnetic NFL near F-QCP to a field-induced F state, 
and we established a magnetic phase diagram of the (Fei_ x Co x )3Mo3N system, where the F 
order is stabilized by a slight doping of the nonmagnetic element. Neutron scattering mea- 
surements suggested that the AF and the F correlations coexist in FesMosN. We proposed 
that the SQ lattice is a new geometrically frustrated system, in which the F and the AF 
interactions compete to select frustrated and non-frustrated states. We applied this model 
to explain the suppression of long-range order in pure FesMosN. In other words, the emer- 
gence of the F-QCP is interpreted as being the result of frustration. Although a number 
of 77-carbide-type transition- metal compounds are known to exist, their quantum physical 
properties have been less studied. We believe that these compounds are good candidates for 



testing geometric frustration in metallic magnets and searching for new and exotic quantum 
phenomena. 
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FIG. 1: (color online) Schematic view of the stella quadrangla lattice. The unit, a stella quadran- 
gula, consists of two nested regular tetrahedra with the same center of gravity. 
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FIG. 2: (color online) (a) High-field magnetization curve of FeaMosN at 4.2 K. The inset shows 
dM/dH in both field- increasing and decreasing processes. The top half of the inset shows a 
magnified version of the image in the bottom half, (b) T dependences of the metamagnetic field 
Hq (average value of field-increasing and decreasing processes, closed circles) and the field of the 
small anomaly in the field-increasing process, He (open circles), (c) T dependence of the hysteresis 
width AH. 
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FIG. 3: (color online) (a) T dependences of inverse susceptibility and spontaneous magnetization 
(estimated from the Arrott plot) of (Feo.sCoo.2)3Mo3N. (b) Arrott plots of (Feo.sCoo.2)3Mo3N. (c) 
Magnetic phase diagram of the (Fei_ x Co x )3Mo3N system. Tq (open circles) values were determined 
from Arrott plots. Tqm (closed circle) is the hypothetical Tc for pure Fe3Mo3N. 
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FIG. 4: (color online) Wave number {K) dependences of quasi-elastic neutron scattering of 
FesMosN at 5.5 K for energy transfers of E = 1 and 2 meV. K is in the reciprocal lattice unit 
(r.l.u.). The solid curve represents the scattering function S(K,E) calculated for the SQ lattice 
(see text). The broken curve represents the K dependence of the squared magnetic form factor for 
Fe 3+ . The inset shows the E scans measured at K = 1.22, 1.73, and 2.83 in r.l.u. Solid curves 
represent the best fits to the Lorentzian quasi-elastic scattering function. 
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FIG. 5: Typical dispersion curves of J{q) calculated for the SQ lattice with negative J2. (a) J<i- 
dominant i\J\\/ J2 = —0.5) and (b) Ji-dominant (—0.8) cases. The results do not depend on the 
sign of J\. In these calculations, z = 0.2937 was used as the coordinate of 32e, and the boundary 
of these states appears between \J\\j Ji = —0.6 and —0.8. 



13 



